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Abstract. The Majority Rule is applied to a topology that consists of two coupled 
random networks, thereby mimicking the modular structure observed in social 
networks. We calculate analytically the asymptotic behaviour of the model and derive a 
phase diagram that depends on the frequency of random opinion flips and on the inter- 
connectivity between the two communities. It is shown that three regimes may take 
place: a disordered regime, where no collective phenomena takes place; a symmetric 
regime, where the nodes in both communities reach the same average opinion; an 
asymmetric regime, where the nodes in each community reach an opposite average 
opinion. The transition from the asymmetric regime to the symmetric regime is shown 
to be discontinuous. 
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1. Introduction 

In the last few years, the study of networks has received an enormous amount of attention 
from the scientific community [H[2], in disciplines as diverse as computer and information 
sciences (the Internet and the World Wide Web), sociology and epidemiology (networks 
of personal or social contacts), biology (metabolic and protein interaction), etc. This 
outburst of interest has been driven mainly by the possibility to use networks in order 
to represent many complex systems and by the availability of communication networks 
and computers that allow us to gather and analyze data on a scale far larger than 
previously possible. The resulting change of scale (from a few dozen of nodes in 
earlier works to several thousands of nodes today) has not only lead to the definition 
of new statistical quantities in order to describe large networks, e.g. degree distribution 
or clustering coefficient, but it has also addressed problems pertaining to Statistical 
Physics, for instance by looking at the interplay between the microscopic interactions 
of neighbouring nodes and the behaviour of the system at a macroscopic level. Such 
a problem takes place in social networks, i.e. nodes represent individuals and links 
between nodes represent their relations (e.g. friendship, co-authorship), when one tries 
to find macroscopic equations for the evolution of "society". Indeed, many studies have 
revealed non-trivial structures in social networks, such as fat-tailed degree distributions 
[3], a high clustering coefficient [I] and the presence of communities [5]. A primordial 
problem is therefore to understand how this underlying topology influences the way 
the interactions between individuals (physical contact, discussions) may (or not) lead 
to collective phenomena. Typical examples would be the propagation of a virus [6] or 
opinion [TJ [8] in a social network, that may lead to the outbreak of an epidemics or of 
a new trend/fashion. 

It is now well-known that degree heterogeneity [HI [9] is an important factor that 
may radically alter the macroscopic behaviour of a network but, surprisingly, the role 
played by its modular structure is still poorly known [TUJ [11]. It has been observed, 
though, that many social networks exhibit modular structures [HI EEJ 03] > i- e - they 
are composed of highly connected communities, while nodes in different communities 
are sparsely connected. This lack of interaction between communities certainly has 
consequences on the way information diffuses through the network, for instance, but it 
also suggests that nodes belonging to different communities may behave in a radically 
different way. 

In this paper, we address such a problem by focusing on a simple model for networks 
with two communities, the Coupled Random Networks (CRN). To do so, one consider a 
set of N nodes that one divides into two classes and one randomly assigns links between 
the nodes. Moreover, one assumes that the probability for two nodes to be linked is 
larger when they belong to the same class. Let us stress that CRN has been first 
introduced in [5] and that it is particularly suitable in order to highlight the role of the 
network modularity while preserving its randomness. The microscopic dynamics that 
we apply on CRN is the well-known Majority Rule (MR) p3]. MR is a very general 
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model for opinion-formation, i.e. nodes copy the behaviour of their neighbour, thereby 
suggesting that the results derived in this paper should also apply to other models of 
the same family. The effect of the inter-connectivity v and of the frequency of random 
flips, measured by the parameter q (~ temperature) on the phase diagram is studied 
analytically. It is shown that three regimes may take place, depending on the parameters 
and on the initial conditions: a disordered regime, where no collective phenomena takes 
place; a symmetric regime, where the nodes in both communities reach the same average 
opinion; an asymmetric regime, where the nodes in each community reach an opposite 
average opinion. The transition from the asymmetric regime to the symmetric regime 
is shown to be discontinuous. It is remarkable to note that a similar discontinuous 
transition also takes place when one applies MR to another network with communities, 
namely the Coupled Fully- Connected Networks introduced in [TUj. The main advantage 
of CRN is that its simpler structure allows to perform all the calculations exactly and 
to consider the case of a non-vanishing q in detail. 

2. Majority Rule 

The network is composed of N nodes, each of them endowed with an opinion that can be 
either a or (3. At each time step, one of the nodes is randomly selected and two processes 
may take place. With probability q, the selected node randomly picks an opinion a or (3, 
whatever its previous opinion or the opinion of its neighbours. With probability 1 — q, 
two neighbouring nodes of the selected node are also selected and the three agents in 
this majority triplet all adopt the state of the local majority (see Fig.l). The parameter 
q therefore measures the competition between individual choices, that have a tendency 
to randomize the opinions in the system, and neighbouring interactions, that tend to 
homogenize the opinion of agents. In the case q — 0, it is well-known that the system 
asymptotically reaches global consensus where all nodes share the same opinion [14J. In 
the other limiting case q = 1, the system is purely random and the average (over the 
realizations of the random process) number of nodes with opinion a at time t, denoted 
by A t , goes to N/2. 

Let us first focus on a network of individuals that are highly connected (in order to 
justify the use of mean-field methods) and where all the nodes are equivalent. That case 
has been studied in detail elsewhere [T5] . and is repeated here for the sake of clarity and 
introducing notations. It is straightforward to show that the mean-field rate equation 
for A t reads 



where a t = A t /N is the average proportion of nodes with opinion a. The term 
proportional to q accounts for the individual random flips. The second term, accounting 
for majority processes, is evaluated by calculating the probability that the majority 
triplet is composed of two nodes a and one node (3, 3af(l — a t ), or of two nodes (3 and 




(1) 
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Figure 1. Sketch of a time step, where one node (surrounded in black) and two of its 
neighbours (surrounded in grey) are selected. The majority rule implies that the blue 
node become red. 



one node a, 3a t (l — a t ) 2 . Consequently, the total contribution to the evolution of A t is 

3 (a 2 t (l - a*) - a t (l - a t f) = -3a t (l - 3a t + 2a 2 t ). (2) 

It is easy to show that a = 1/2 is always a stationary solution of Eq.([T]), as expected 
from symmetry reasons, a = 1/2 corresponds to a disordered state where no collective 
opinion has emerged in the system. It is useful to rewrite the evolution equation for the 
quantities A t = A t - N/2 and 8 t = A t /N = a - 1/2 



A 



5; 



t+1 



A t + j(3-5q- 12(1 - q)5 2 t 



(3) 



from which one finds that the disordered solution a = 1/2 ceases to be stable when 
q < 3/5. In that case, the system reaches one of the following asymmetric solutions 
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The system therefore undergoes an order-disorder transition at q = 3/5. Under this 
value, a collective opinion has emerged due to the imitation between neighbouring nodes. 
In the limit case q — > 0, one finds a_ = and a + = 1 in agreement with the results of 

Em. 



3. Role of communities 

3.1. Coupled Random Networks 

Distinct communities within networks are defined as subsets of nodes which are more 
densely linked when compared to the rest of the network. For the sake of simplicity, 
we restrict the scope to networks composed of only two communities, denoted by 1 and 
2. Our goal is to build an uncorrelated random network where nodes in 1 are more 
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Figure 2. Typical realizations of coupled random networks for small values of v. The 
network is composed of N = 100 nodes and pi„ =0.1. The system obviously consists 
of two communities that are less and less discernible for larger values of v. The graphs 
were plotted thanks to the visone graphical tools [16] . 



likely to be connected with nodes in 1 than with nodes in 2, and vice-versa. To do 
so, we consider a set of N nodes that we divide nodes into two classes, 1 and 2. We 
evaluate each pair of nodes in the system and draw a link between these nodes with 
probability py, where i and j G {1, 2} are the class to which the two nodes belong. In 
the following, we consider a case where the number of nodes in 1 and 2, respectively 
denoted by N\ and N 2 , are equal N\ = N 2 = N/2. Moreover, we restrict the scope 
to the following probabilities p\ 2 = P21 = Pcross and pn = p 22 = Pin- By construction, 
nodes in 1 are therefore connected on average to ki n = Pi n (N — l)/2 ~ pi n N/2 nodes in 
1 and to k cross = p cross N/2 nodes in 2, while nodes in 2 are connected to k cross nodes in 
1 and k in nodes in 2. Let us stress that this approach pertains to the theory of networks 
with hidden variables [TF1 [18] where the probability for a pair to be linked would not 
factorize Pij 7^ PiPj. A similar model has already been used in order to test methods for 
community detection in [5]. 

This binary mixture, that we call a Coupled Random Network (CRN) is particularly 
suitable in order to reveal the role of network modularity Indeed, the inter- 

connectivity between the communities is tunable through the parameter v = p cr0 ss/ ' Pin- 
In the following, we focus on the interval v e [0, 1] for which inter-community links are 
less frequent than intra-community links. When v — *■ 1, one recovers an homogeneous 
network where all nodes are a priori equivalent, while when v « 1, the communities 
are sparsely connected with each other. 

Before going further, one should also point to an alternative model of modular 
networks introduced in [10]. In that paper, Coupled Fully Connected Networks are 
composed of three kinds of nodes, the interface nodes, the nodes of type 1 and the 
nodes of type 2, and an inter-connectivity parameter measures the modularity of the 
network. MR has been applied on this topology with q = and it has been shown that 
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a discontinuous transition from an asymmetric to a symmetric state takes place. In the 
following, we will not only show that a similar discontinuous transition takes place on 
CRN, but we will also study analytically the behaviour of the system for q ^ 



3.2. Equation of evolution 

Let us denote by A\ and A 2 the average number of nodes with opinion a among the 
two types of nodes. Let us first focus on the contributions when majority triplets are 
selected. At each time step, the probability that the selected node belongs to the first 
community is 1/2. In that case, the probability that a randomly chosen link around 
the selected node goes to a node in 1 is k in /(k in + k cross ) = 1/(1 + u). The probability 
that this randomly chosen link goes to a node in 2 is k cross /(ki n + k cross ) = u/(l + v). 
Consequently, the probability that the selected node belongs to 1 and that both of its 
selected neighbours belong to 2 is 

1 t (5, 
2(1 + ■/)*• (5) 

Similarly, the probability that the selected node belongs to 1, that one of its neighbours 

belongs to 1 and that the other neighbour belongs to 2 is 

1 2v 



2(l + z/) 2 ' 

while the probability that all three nodes belong to 1 is 
1 1 



(6) 



(7) 



2(1 + v) 2 ' 

The probabilities of events when the selected node belongs to 2 are found in the same 
way. Putting all contributions together, one finds the probabilities P( x , y ) that x nodes 
1 and y nodes 2 belong to the majority triplet 

1 1 
2(l + z/) 2 

1 2u 1 v 2 1 v 2 + 2v 



P (3,0) 



P(2,l) = - ~ . ^ + - 

^(1,2) 



2(l + z/) 2 2(1 + z^) 2 2(l + i/) 
1 v 2 + 2u 



2(l + u) 2 

P(0 < 3) = 2 oTW (8) 

where the normalization J2 xy P{x,y) = 1 is verified. In order to derive coupled equations 
for Ai- t and A 2 -t that would generalize Eq.flT]), one needs to evaluate the evolution of 
these quantities when a triplet (x, y) is selected. To do so, one follows the steps described 
in [10] and, when q — 0, one obtains the equation of evolution 



A 1;t+1 - Ax-t = ^^ 2 ( Q i & i ~ a i & i) + Jl ^p ( fl 2fli&i - 016261) 

J- v ~T / 2r 7 2\ 
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2 (l + vy (l + vy 

1 v 2 + 2v , 9 , , 9s 
+ 2 (l + z/)2 (a ' &2 ~ a2 ^' (9) 

where Oj and 6j are respectively the proportion of nodes with opinion a and /? in the 
community i. After incorporating the term due to random flips, proportional to q, one 
obtains the set of non-linear equations 

A l]t+l - A l;t = I - + (1 - ^ JT+^f ^ 1 ~ aib ^ 

+ (1 + u ^ a2aibl ~ + 2 (T+^F ~~ 2 ^ 

^ 2 ;i+i - A 2 . t = - - ^ + (1 - — ^ (agfea - a 2& 2 ' 



4 2 ^ H,[ 2(l + v) 2 ^ 2U2 " 2 ^ 

i/ 2 + 2v . , , , . 1 z/ 2 + 2z/ . 2l , . 

+ 7— — r^(aia 2 o 2 - a 2 b x b 2 ) + -7— — rn-( a i & 2 - e^oJJ. (10) 



Direct calculations show that Eq. lflOl) reduces to Eq.fllj) in the limit v = 1, as expected 
due to the indistinguishability of the nodes in that case. 

3. 3. Stability of the disordered solution 

It is straightforward to show that a± = 1/2, a 2 = 1/2 is always a stationary solution 
of Eq. liTOl . whatever the values of v and q. This solution consists of a disordered state 
where both communities behave similarly and where no favourite opinion has emerged 
due to MR. We study the stability of this disordered state by looking at small deviations 
£1 = a\ — 1/2 and e 2 = a 2 — 1/2 and keeping only linear corrections. In the continuous 
time limit and after re-scaling the time units, the evolution equations for these deviations 
read 

-3 + 5g + 2z/(l + g) + z/ 2 (l + g) z/(2 + v)(l - q) 

dtei ~ W+W ei + (1 + vf 62 

v{2 + v){l-q) -3 + 5q + 2is(l + q) + v 2 (l + q) 

dt62 = (i + v) 2 61 4(TT^ e2 ' 

The eigenvalues of this linearized matrix of evolution are 
(3 - 10z/ - 5z/ 2 - 5q + 6uq + 3u 2 q) 
1 ~ 4(1 + z/) 2 

A 2 = i(3-5g). (12) 

By definition, the disordered solution is stable only when both eigenvalues are negative 
[j~9] . thereby ensuring that a small perturbation asymptotically vanishes. It is easy to 
show that only the values of q in the interval ]3/5, 1] respect this condition, whatever 
the value of v. This implies that the location of the order-disorder transition is not 
affected by the modularity of the network. Let us also stress that that \\ goes to A 2 



Coexistence of opposite opinions in a network with communities. 



8 



wmwtowwiww: — I 









W ■ 






z/ = 0.05 














Asymmetric 
Symmetric K 



0.2 



0.4 0.6 

q 



0.8 



Figure 3. Computer simulations of MR on Coupled Random Networks with N = 10 4 
nodes and pin — 0.01. The simulations are stopped after 10 3 steps/node and the results 
averaged over 100 realizations of the random process. The vertical dashed lines point to 
the theoretical transition point obtained from Eq. ([25|) and to the critical value q = 3/5 
(Fig. 3d). The solid lines correspond to the theoretical predictions (fT6|) and (fl~8|) . The 
simulations are either started with a symmetric initial condition a\ = 1, ai = 1 or with 
an asymmetric initial condition a\ = 1, a% = 0. (a) Bifurcation diagram of \a\ — a,2\ 
as a function of for simulations starting from asymmetric initial conditions. The 
system ceases to be asymmetric \a\ — a,2\ > above v c « 0.15. (b) and (c) Bifurcation 
diagram of a\ as a function of v, starting the simulations from asymmetric or symmetric 
initial conditions for q = (b) and q — 0.2 (c). (d) Bifurcation diagram of a\ as a 
function of q, starting the simulations from asymmetric or symmetric initial conditions 
for v — 0.05. In that case, the system may undergo two transitions: one from the 
asymmetric to the symmetric state at q ~ 0.485, and one from the symmetric to the 
disordered state at q = 3/5. 



when v = 0. This is expected, as the system is composed of two independent networks 
in that case, so that the equations of evolution for A\ and A 2 are decoupled. 

3.4- Symmetric solution 

Our knowledge of the case v = 1 (section 2) and the natural symmetry of CRN suggests 
to look at solutions of the form a\ = 1/2 + 5s, a 2 = 1/2 + 5s (S for symmetric). By 
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inserting this ansatz into Eq.( 



q + 8s) ?)[ _3 (§ + ~ ^s) - (| - + fc) 



4 2 ' v " " /L 2 (l + u) 2 

+ ,-, I 3o ((o + 5 s) (o - 5 s) " (tt " 5s) 2 (7T + 5 s)) 



direct calculations lead to the relation 

_ ^ + { k - 251) ^- q h- + + M = 0- (14) 

Let us insist on the fact that Eq.( fl4l) is exact and not an expansion for small 5s- It 
is a direct consequence of the above symmetry ansatz. The disordered solution 5s = 
obviously satisfies Eq.(fl4l). but symmetric solutions are also possible if they satisfy 



?(1 + ^) 2 . ,1 OX2U 3 . 3 .2 



+ (o-2^)[- + -^ + 3H = 0, (15) 



2(1 -q) v 2 bn 2 2 
so that symmetric solutions have the form 

Si = 3 ~ 59 - (16) 
s 12(1 -g) V ; 

It is remarkable to note that the symmetric solution does not depend on the inter- 
connectivity v. This is checked by comparing (Tl6|) with the solution (j3J) obtained when 
the system is composed of only one community. It is also straightforward to show that 
( TT6l) is stable when q < 3/5, as expected, from which one concludes that none of the 
characteristics of the order- disorder transition have been altered by the modularity of 
the network. 

3.5. Asymmetric solution 

Let us first focus on the case q = where the dynamics is only driven by MR. We have 
shown above that the system may reach a symmetric frozen state ai — 1, a<i = 1 or 
Oi = 0, a>2 = in that case (see Eq. lTIBl) ). However, computer simulations (Fig.3a) show 
that an asymmetric stationary state may prevail for small enough values of v. Computer 
simulations also show that the asymmetric state is characterized by averages of the form 
a\ = 1/2 + 5a and a 2 = 1/2 — 5a (A for asymmetric). Based on these numerical results, 
we look for stationary solutions of Eq.([9]) having this form. It is straightforward to show 
that the equations for A\ and Ai lead to the following condition 

(35 A - 125 \) + (v 2 + 2v)(-26 A + 85 3 A ) - (v 2 + 2v)(35 A + 4<%) = 0, (17) 
whose solutions are either 5a = (disordered state) or 
2 _ 3 - 5(^ 2 + 2v) 

§A ~ l2-A(v 2 + 2vy (18) 
The disordered solution 5a = is unstable when q = 0, as shown above, and it it 
thus discarded. Let us also stress that the asymmetric solution differs from the frozen 
symmetric solution by the fact that fluctuations continue to take place in it. 
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By construction, the asymmetric solution exists only when 5\ > 0, namely when 
v E [0, -s+v^o m 0.26]. In order to check the stability of ffl8|) in this interval, we focus 
on the small deviations ei = a\ — (1/2 + 5a) and e<i = a<i — (1/2 — 5a)- After inserting 
these expressions into Eq. (Q and keeping only linear terms, lengthy calculations lead to 
the eigenvalues 

3 - lOv - 5z/ 2 

1_ 2(1 + i/) 2 

A - 3 1 - 6 "- 3 " 2 (1 g) 

from which one shows that the asymmetric solution loses its stability at a critical value 

Uc = (V48-6)/6 « 0.15. (20) 

Consequently, the system exhibits a discontinuous transition at z/ c , as 1/2 + 5a(v c ) 
0.85 7^ 1. When v < i> c , the system may reach either the symmetric or the asymmetric 
solution depending on the initial conditions (and on the fluctuations). When v > v c , in 
contrast, only the symmetric solution is attained in the long time limit. Let us stress that 
MR also undergoes a discontinuous transition from an asymmetric state to a symmetric 
state when it is applied to Coupled Fully Connected Networks [10] . This similarity 
suggests therefore that such a transition is a generic feature of networks with modular 
structure. The above theoretical results have been verified by performing computer 
simulations of MR. The asymmetric solution (1181) and the location of the transition ( 1201) 
are in perfect agreement with the simulations (Fig. 3b). 

It is straightforward to generalize these results when q > 0. To do so, one inserts 
ai = 1/2 + 5a, CL2 — 1/2 — 5a into Eq. tfroj) from which one obtains the relation 



- 2q5 A + q l m A - 125\) + (y 1 + 2u)(-25 A + 85 3 A ) 
(1 + vy 

-(u 2 + 2u)(35 A + 4:5 3 A )} = 0. (21) 
The stationary solutions of (I2T1) are either 5a = or 



2 3 _ 2^gL _ 5(^2 + 2z/) 
6a = 12 -4(1/2 + 2v) > ( 22 ) 
and the eigenvalues of the linearized equation of evolution around the stationary solution 
are 

3 - 10i/ - 5z/ 2 - 5q + 6uq + 3u 2 q 

1_ 2(1 + z/) 2 

3 - 18u - 9u 2 -5q+ lAvq + 7u 2 q 

2(l + i/) 

The asymmetric solution is shown to loose its stability when 

3 — 18i/ — 9z/ 2 - 5q + 14i/g + 7z/ 2 g = (24) 
that one simplifies into 

^(g) = -l + 2ySH| (25) 



A 2 0/1 , .x 2 • 
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Figure 4. Phase diagram of MR on CRN. Three phases may take place, i) a disordered 
phase when q > 3/5; ii) a symmetric phase when q < 3/5; iii) an asymmetric phase 
when q < 3/5 and when v < — 1 + 2>J j^Zg ■ A system in the asymmetric state, but 
close to the transition line, e.g. the green triangle, may lose its stability due to an 
increase of the number of inter-community links (along v) or to an increase of the 
internal fluctuations (along q). 



This relation therefore determines the critical line above which only symmetric solutions 
prevail (see Fig.4). One can shown that (1251) decreases with q and that it goes to zero 
at the transition point q = 3/5. It is also easy to show that the transition from the 
asymmetric to symmetric state is discontinuous for any values of q < 3/5. Eq. f l22l) 
and Eq. (125]) have been successfully checked by computer simulations (see Fig. 3c and 
Fig.3d). 



4. Discussions 



In this last section, we would like to point to some possible applications of this work and 
to the implications of our theoretical predictions. First of all, let us remind that many 
networks exhibit modular structures. It is therefore of significant practical importance 
to better understand the role played by these structures. This is true for social networks, 
where groups of nodes correspond to social communities or cliques, but also for the World 
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Wide Web, where groups of highly connected sites may be related by a common topic 
[22], [23] . biological networks, where clusters of nodes may be associated to functional 
modules [2U [25] [26] , GDP networks, from which dependencies between countries might 
be uncovered [271 [28] and even self-citation networks where clusters may reveal the 
field mobility of a scientist [29]. The identification of such communities is thus a 
very important problem, that has received a considerable amount of attention recently 
[30] I3T] 1321 [33] 134"] [35] [36] [37] . Amongst others investigations, it has been shown that 
the use of Ising-like models may be helpful in order to unravel communities in complex 
networks [38] [39] • This approach consists in applying a dynamics such as the majority 
rule on a pre-given network and to identify a community as a set of nodes having the 
same opinion/spin. In the language of this paper, the identification of the communities 
is possible only if the system reaches an asymmetric regime, i.e. if both communities 
reach different states. Our work could therefore provide a theoretical background for 
the use of such identification techniques. For instance, our results show that MR has 
a "minimum resolution", i.e. MR does not discriminate communities when v > v c (q) 
(with u c (q) < (V^ — 6)/6) because the asymmetric state is not stable in that case. One 
expects that similar effects could also take place for similar models of opinion formation. 

The fact that nodes in different communities exhibit different behaviours, i.e. the 
asymmetric regime, has been observed in many situations. Amongst others, one may 
think of sub communities in collaboration networks, that may correspond roughly to 
topics of research [5], the Political Blogosphere, where it was observed that bloggers 
having a different political opinion are segregated [10], the existence of niche markets 
[41] . where small communities may use products that are different from those used by 
the majority, language dynamics, where it is well-known that natural frontiers may 
also coincide with a linguistic frontier [12], etc. The discontinuity of the transition 
from this asymmetric state to a symmetric "globalized" state might thus have radical 
consequences for systems close to the transition line: the addition of a few links between 
the communities or a small increase of the fluctuations inside the system (see Fig.4) may 
be sufficient in order to drive the system out of the asymmetric state. Such rapid shifts 
of behaviour of a whole sub-community should be searched in empirical data, in the 
dynamics of trends or fashion [13] for instance. 

To conclude, we would like to point to a possible generalization that would certainly 
make the model more realistic. In this paper, we have focused on a model of opinion 
formation evolving on a static topology. The effect of the network structure on the 
formation of a collective opinion has therefore been highlighted, but the fact that 
the opinion of the nodes themselves might influence the topology around them, e.g. 
links between nodes with different opinions could disappear, has not been considered. 
A model where the network inter-connectivity v might co-evolve |44j with the node 
spins/opinions would thus be of high interest and will be considered in a future work. 
Other interesting generalizations would incorporate agents with ageing [l5l H6] H7J or 
memory [IB], the presence of leaders [10], etc. 
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